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ABSTRACT 

We examine high-resolution gravitational N-body simulations of the fl = 1 cold dark 
matter (CDM) model in order to determine whether there is any normalization of the 
initial density fluctuation spectrum that yields acceptable results for galaxy clustering and 
velocities. Dense dark matter halos in the evolved mass distribution are identified with 
luminous galaxies; the most massive halos are also considered as sites for galaxy groups, 
with a range of possibilities explored for the group mass to light ratios. We verify the 
earlier conclusions of White et al. (1987) for the low amplitude (high bias) CDM model 
— the galaxy correlation function is marginally acceptable but that there are too many 
galaxies. We also show that the peak biasing method does not accurately reproduce the 
results obtained using dense halos identified in the simulations themselves. The COBE 
anisotropy implies a higher normalization, resulting in problems with excessive pairwise 
galaxy velocity dispersion unless a strong velocity bias is present. Although we confirm 
the strong velocity bias of halos reported by Couchman & Carlberg (1992), we show that 
the galaxy motions are still too large on small scales. We find no amplitude for which 
the CDM model can reconcile simultaneously the galaxy correlation function, the low 
pairwise velocity dispersion, and the richness distribution of groups and clusters. With 
the normalization implied by COBE, the CDM spectrum has too much power on small 
scales if O = 1. 



Subject headings: cosmology: theory — dark matter — galaxies: clustering — galaxies: 
formation 



1. INTRODUCTION 



The cold dark matter (CDM) model of galaxy formation has had a checkered history. 
First proposed by Peebles in 1982, the model has the virtue of being relatively well- 
defined and testable. Assuming Hq = 50kms~^ Mpc~^ and = 1, the only fundamental 
free parameter is the overall amplitude of density fluctuations. We characterize this 
by the conventional quantity cTg, defined to be the rms density fluctuation, using the 
linear power spectrum, in a sphere of radius 800kms~^. (See, e.g., Bertschinger 1992 
for discussion of this and alternative conventions for the normalization.) Once erg is 
specified the CDM model has, in principle, strong predictive power, although many of 
the predictions require N-body and dissipative numerical computations. The complexity 
of the nonlinear evolution and dissipation has led over the last decade to a lively debate 
concerning the viability of the CDM model. 

In 1985, Davis et al. showed that the CDM model cannot simultaneously fit galaxy 
clustering (i.e., two-point correlation function) and small-scale velocities (i.e., pairwise 
velocity dispersion) for any erg if dark matter traces galaxies. For erg = 1, motivated by 
the observation that ug ~ 1 for galaxies (Davis & Peebles 1983), the relative velocities 
of galaxies are predicted to be much larger than observed. The solution proposed by 
Davis et al. was to decrease the amplitude of density fluctuations by a factor 2.5 to 
(Tg = 0.4, thereby decreasing the pairwise velocity dispersion of galaxies to roughly match 
observations. In the process, the clustering of galaxies was also diminished. The two- 
point correlation function was boosted up to the observed range by assuming that galaxies 
form only in the initially highest-density regions, according to the "peak biasing" scheme 
proposed by Kaiser (1984). Roughly speaking, galaxy density fluctuations were assumed 
to be 2.5 times the dark matter fluctuations (although peak biasing does not give an 
exactly linear relation between galaxies and mass fluctuations). The stronger correlations 
introduced with biasing compensated for the smaller correlations (and velocities) resulting 
from lowering the amplitude to erg = 0.4. The paper of Davis et al. (1985, hereafter 
DEFW) marked the birth of the biased CDM model. 

The amplitude of the CDM density fluctuations affects large-scale structure as well 
as small-scale (20 Mpc and less) clustering and velocities. Many authors have pointed out 
that the low amplitude of biased CDM is in apparent conflict with large-scale structure. 
The freedom to vary the biasing (ratio of galaxies to mass) makes it somewhat difficult to 
pin down these problems. However, large-scale peculiar velocities arc particularly useful 
for testing the normalization because galaxies should trace the same large-scale flows as 
does dark matter — all bodies accelerate the same way in a gravitational field. Motivated 
by this fact, Bertschinger & Juszkiewicz (1988) tested CDM predictions against the "great 
attractor" fits of Faber & Burstein (1988) and concluded that CDM with erg < 2/3 was 
inconsistent with the data. This conclusion was strengthened by analysis of several large- 
scale galaxy surveys: the angular correlation function of Maddox et al. (1990) from the 
APM survey; the moments of galaxy counts in cells of Saunders et al. (1991) from the 
IRAS/QDOT survey; and the galaxy density power spectrum from the CfA2 redshift 
survey by Vogeley et al. (1992). 

The large-scale structure measurements indicate a need for more power, hence a 
larger ug if the CDM power spectrum is retained. However, high- amplitude CDM faces 
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the problem of large velocity dispersion noted originally by Davis et al. (1985). A possible 
solution was suggested by Carlberg & Couchman (1989) and Carlberg, Couchman, and 
Thomas (1990): velocity bias. They noted that the pairwise velocity dispersion of dark 
matter halos in high-resolution N-body simulations is substantially less than that of the 
mass. This effect was not discovered by Davis et al. because their simulations did not 
have the resolution needed to find galaxy halos composed of many particles (although it 
could have been found by White et al. 1987). Couchman & Carlberg (1992) pointed out 
that with an amplitude corresponding to as = 1, CDM would do well on large scales, 
while clustering and velocity biasing might solve the problems on small scales. 

However, it is not enough for CDM to predict the correct two-point correlation 
function and pairwise velocity dispersion of galaxies. It must also predict the correct 
abundance of galaxies and of galaxy groups as a function of richness. Testing these 
requires high-resolution numerical simulations. An important first step was taken in 
1987 by White et al. (hereafter WDEF). They performed a P^M simulation, evolved to 
as — 0.4, in a 50 Mpc box with enough particles (64^) to study resolved dark matter 
halos. They found that the evolved halos are, indeed, more strongly correlated than the 
mass, with the correlation function in reasonable agreement with observations for halos 
with circular rotation speeds exceeding 250kms~^. However, they found the numbers of 
halos in their simulations with Vch-c ^ lOOkms"^, after breaking up the overly-merged 
massive halos, to be greater than the observed numbers by a factor of 2 or more. 

The uncertainty over the correct amplitude for normalizing the CDM (or any other) 
power spectrum has largely ended with the measurement of the cosmic microwave back- 
ground anisotropy by Smoot et al. (1992). Their measurements imply ag ~ 1.1 (Wright 
et al. 1992; Efstathiou, Bond, & White 1992; Adams et al. 1993) if the spectrum is 
scale-invariant as we suppose. Consequently, all standard CDM models with dg — 0.4 
are obsolete. On the other hand, CDM with as = 1 might be an attractive model if the 
small-scale velocity bias found by Carlberg et al. is sufficiently strong and if the galaxy 
numbers and group multiplicities can be made reasonably to match the observations. 

In previous work (Bertschinger & Gelb 1991) we have also found a strong velocity 
bias for halos in the CDM model. However, our interpretation differs somewhat from 
that of Carlberg (1991); the reduction appears to be a statistical effect, arising because 
the pairwise velocity statistic weights galaxies by pairs and quadratically by velocities. 
Therefore, when it is applied to all mass particles, the massive halos in galaxy clusters, 
with velocities comparable to the cluster dispersion, contribute strongly to the overall 
pairwise velocity dispersion. When the simulated halos are used, on the other hand, 
strong merging in the clusters eliminates most of the halos and leaves a single massive 
object in the center, which has little weight in the pairwise sum. The small number of 
halos in clusters reduces the pairwise velocity dispersion of halos, but it also leads to 
clusters with far too few objects to be identified as galaxies. 

In a preceding paper (Gelb & Bertschinger 1993, hereafter Paper I), we explored 
in detail the distribution of simulated halos by circular velocities, and we concluded 
that no value of erg could fully satisfy the constraints given by the number density of 
galaxies. For any value of ug one had too many halos, assuming that circular velocity can 
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be related to luminosity by TuUy-Fisher and Faber-Jackson relations. However, these 
results are limited by the fact that the N-body simulations we used lack gas, so it is 
worthwhile making other tests of the model, with generous allowance for uncertainties in 
how galaxies are related to dark matter halos. 

Dissipationless simulations like ours have the disadvantage of being unable to cor- 
rectly model the process of galaxy formation inside dark matter halos. By necessity, we 
assume that galaxies form only inside dark matter halos. This assumption appears to 
be confirmed by gas dynamical simulations (e.g., Katz, Hernquist, and Weinberg 1992; 
Cen & Ostriker 1992a,b; Evrard, Summers, &; Davis 1994), which also support our con- 
clusions concerning the degree of velocity biasing. Why then should we continue with 
dissipationless simulations? The reason is dynamic range. With simulations including 
gravity only, we are able to resolve galaxy halos in volumes up to 100 Mpc on a side, large 
enough to capture the long wavelength density fluctuations important for high- amplitude 
CDM (cf. Paper I). Gas dynamical simulations with equal numbers of particles are still 
prohibitively expensive. The volumes studied to date with high resolution gas dynamical 
simulations have been too small to include all important dynamical effects. 

In this paper we study spatial clustering and velocity statistics of dark matter halos in 
the CDM model using high-resolution N-body simulations. The principle goal is to answer 
the question: is there a normalization erg such that the two-point spatial correlation 
function of the resolved halos and the pairwise velocity dispersion of the resolved halos 
matches the observations? This question is addressed by analyzing the particle-mesh 
and particle-particle/particle- mesh (P^M) N-body simulations discussed in paper I. For 
economy of notation (see Paper I) we refer to the simulations as CDMn(A^'^,L,J?i/2)- The 
numbers in parentheses indicate the following simulation parameters: 1) particles, 2) 
a comoving box of length L Mpc on a side, and 3) a comoving force softening length of 
Ri/2 kpc. The force softening is characterized by r = Ri/2'- where r'^Fr/{Gm?) = 1/2, 
i.e. half its Newtonian value. 

Our studies focus on the simulation CDM16(144^,100,85). The comoving Plummer 
softening length is e = 65 kpc and the particle mass is mpart = 2.3 x IO^^Mq. Of all our 
simulations (see Gelb 1992, Appendix II), this simulation has the best compromise of mass 
and force resolution in a 100 Mpc box. A relatively large box is required to adequately 
represent waves in the initial conditions, particularly for evolution up to ds = 1 (see 
Paper I and §§ 2.1 and 2.2 below). We will examine positions and peculiar velocities 
at several amplitudes: as — 0.4, 0.5, 0.7, and 1.0. Each epoch studied is considered a 
candidate for the present day; i.e. we test whether spatial and velocity statistics match 
the observations. The box length is assumed to have a physical length of 100 Mpc for 
CDM16 with a Hubble constant Hq = 50kms~-^ Mpc~^ (and f2 = 1) at each candidate 
epoch. 

We will begin by exploring some background material: box size, the standard CDM 
model, massive halos, and velocities in § 2. We then briefly discuss some limitations of 
the method of peak particles and argue the necessity of using resolved halos to study 
the CDM model in § 3. In § 4 we study the two-point correlation function of simulated 
galaxies, while in § 5 we study the abundance and richness distribution of galaxy groups. 
In both cases the results depend on how over-merged halos are broken apart, a point that 
we investigate in some detail. In § 6 we investigate the small-scale velocity statistics of 
the galaxies. In § 7 we summarize the implications for the = 1 CDM model. Further 
numerical details can be found in Bertschinger & Gelb (1991), Gelb (1992), and Paper I. 
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2. BACKGROUND AND THE STANDARD MODEL 



In this section we explore the issue of box size; we test whether we can reproduce the 
results of White et al. (1987); and we consider complications arising from massive halos 
and from definitions of halo velocities. We will confirm that the "standard" biased CDM 
model {Q — 1, as = 0.4) produces far too many halos compared with the observations. 

2.1. Box Size: Mass Correlation Length 

We are interested in studying models evolved to erg = 1. More highly evolved models 
require larger boxes since successively larger waves begin to go nonlinear. In this subsec- 
tion we examine linear theory predictions regarding the importance of long waves in the 
initial conditions. We then present nonlinear results using evolved A'^-body simulations to 
study the dependence on box size of the two-point correlation function and the pairwise 
velocity dispersions of the mass. 

In Figure la we show the linear, rms mass fluctuation in a sphere of radius 8h~^ = 
16 Mpc. We normalize the Holtzman (1989) (5% baryons) CDM power spectrum so 
that (Tg = 1. We then compute crg(Ajnax) by including only waves with wavelength less 
than Amax in the numerical integration of erg. The value of the maximum wavelength 
represented in a simulation computed in a box of length L on a side is Amax — L. As 
Amax — oo, crg(Amax) — ^ 1 by definition. We show vertical bars at Amax = 51.2 Mpc and 
100 Mpc, the sizes of several of our simulation cubes. 

For Amax = 51.2 Mpc we find a-g(Amax) ~ 0.6. For Amax = 100 Mpc we find 
C8('^max) ~ 0.9. For Amax = 150 Mpc wc find crg(Amax) ~ 0.98 and it quickly approaches 
unity thereafter. We conclude that long waves (with A > 50 Mpc) make significant con- 
tributions to CTg and, by extension, to the correlation length tq (where the two-point 
correlation function ^ = 1). Therefore we expect that the nonlinear evolution of our 
A'"-body simulations up to the amplitude erg = 1 will result in a serious underestimate of 
ro in 51.2 Mpc boxes but probably not in 100 Mpc boxes. 

To test the above prediction, in Figure 2 we show ^(r) for various simulations for 
the mass (particles). The important parameter for our discussion is the box size. The 
plots show that the correlation length tq grows roughly linearly with ug from 5 Mpc at 
CTg = 0.5 to 10 Mpc at CTg = 1 for the simulations in larger boxes. We must match the 
observations with the simulated halos, not the mass, unless the halos trace the mass. 
However, we expect the same waves that affect the mass will also affect halo clustering. 
We notice in Figure 2 that the 51.2 Mpc box simulations underestimate the correlation 
length (by about 30% for cxg = 1) when compared with the simulations in > 100 Mpc 
boxes. The underestimate is greater for increasing erg (being about 20% for erg = 0.5), 
as expected. Boxes smaller than 50 Mpc on a side are too small to get the correlation 
length correct to better than 10% even for an amplitude as small as erg = 0.4. However, 
we find that a 100 Mpc box is sufficiently large, as the correlation length (10 Mpc) agrees 
with the results from the 400 Mpc box simulation, despite the poor resolution (2 Mpc) 
of the latter. 

From Figure 2 we can also examine the the simulation-to-simulation variation in the 
value of tq. We see basic agreement in ro among the larger box simulations. However, we 
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do see a smaller value in ro at ug = 0.7 for CDM16 (100 Mpc box) indicating that there 
are simulation-to-simulation fluctuations. The fluctuations for the five R1/2 = 280 kpc 
PM simulations in 51.2 Mpc boxes are shown with la error bars. The largest fluctuations 
are found on the largest scales. 

We conclude that the 51.2 Mpc boxes are too small for accurate predictions of the 
two-point correlation function particularly for larger values of erg. It is unfortunate that 
CDM16 (100 Mpc box) has a slight sag on scales near 10 Mpc owing to a statistical fluc- 
tuation. This will limit some of the conclusions we can draw from this single simulation. 
Nevertheless, the large numbers of halos in this large box provide a fair test of the CDM 
model. More realizations would be helpful, but our other ~ 100 Mpc box simulations 
have poor force and mass resolution. 

2.2. Box Size: Mass Pairwise Velocity Dispersion 

In this subsection we examine the effect of the box size on the pairwise velocity- 
dispersion (Tp(r) as a function of galaxy separation: 

<^l{r)^l{{v2-v,)') , (2.1) 

where the average is taken over pairs of particles (for the mass) or halos separated by 
distance r = |rL — r2 |, with peculiar velocities vi and V2. Nonlinear studies of velocity 
dispersions are the subject of Gelb, Gradwohl, & Frieman (1993). 

For an initial estimate of effects of flnite box size wc consider the linear theory 
prediction for the three-dimensional pairwise dispersion (T„(r) = 3^/^o"p(r). We evaluate 
this quantity using the as = 1 linear normalization, with a wavelength cutoff Amax applied 
to the numerical integration (as we did in § 2.1) to mimic the effect of a finite box size. 
The results are shown in Figure lb. We see there is a significant difference between a 25 
Mpc box and a 50 Mpc box (particularly for larger values of r) and a smaller difference 
between a 50 Mpc box and a 100 Mpc box. The results converge on the scales of interest 
in this paper, i.e. r ~ 1 — 10 Mpc, for A^ax ~ 100 Mpc. 

In Figure 3 we show a"p(r) for the mass for the same simulations studied in § 2.1. We 
notice that the velocities are higher for the simulations in boxes of size > 100 Mpc than for 
the simulations in 51.2 Mpc boxes — long waves in the initial conditions affect nonlinear 
pairwise velocities on smaller scales (see also Gelb, Gradwohl, & Frieman 1993). We see 
that 100 Mpc boxes are sufficient because the results for the 400 Mpc box simulation 
are comparable to the other 100 Mpc box simulations. The 400 Mpc box simulation has 
extremely soft forces {Ri/2—2 Mpc) so the velocities are actually lower than in the 100 
Mpc simulations. (The parameters in the 400 Mpc simulation are close to the values used 
by Park 1990.) We use our results from linear theory to strengthen our argument that 
a 100 Mpc box is sufficient for studying the velocities. This is encouraging because we 
demonstrated in the previous section that a 100 Mpc box is sufficient for studying spatial 
correlations. 

2.3. The Standard Biased CDM Model: as = 0.4 
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In the above sections we found that a box of size 51.2 Mpc is too smaU to include aU 
the long-wavelength contributions important for clustering and pairwise velocity disper- 
sions on small scales. However, much of the past work has been done using simulations of 
this size or smaller. To show that our simulations are consistent with previous work, we 
perform two simulations with parameters similar to those used by WDEF. These P^M 
simulations, CDM12 and CDM13, use 64^ particles in 51.2 Mpc boxes with Plummer 
softening e = 40 kpc comoving. The particle mass is mpart = 3.5 x lO^'^ M©. For compar- 
ison, WDEF computed three 64^ particle P"^M simulations in 50 Mpc boxes with force 
resolution ~ 50 kpc comoving (they used a linear sphere density profile). We study the 
models at erg = 0.4, the normalization advocated by DEFW. To identify dark matter ha- 
los we use two different prescriptions: DENMAX (our density maxima finder, see Paper 
I) and EOF (friends-of-friends); WDEF used the latter. For friends-of- friends we report 
the linking length, I, in units of the mean interparticle spacing. We study the numbers 
of simulated halos and we break up the massive halos into clusters to study the effect on 
the spatial clustering of massive halos with circular velocities Vdrc > 250kms~-^. 

In Table 1 we list the numbers of halos with V^irc > 100, 200, and 250 kms~^. 
The results are shown for averages from CDM12(643,51.2,52) and CDM13(643,51.2,52) 
(two different sets of initial random numbers) and for averages from WDEF. (We divide 
their numbers by three since they report totals for three simulations. We also scale their 
numbers in a 50 Mpc box to a 51.2 Mpc box.) In computing observational estimates we 
characterize all our halos by their circular velocities, and we relate observed estimates of 
one-dimensional velocity dispersions, cri, to circular velocities by 

= . (2.2) 



We report our observational estimates for F = 1 and F = 1.1 (see Paper I) and 
those of WDEF for F = 1 (in our notation). WDEF used F = 1 but in a subsequent 
paper at ag = 0.4 they used F = 1.1 (see Frenk et al. 1988). We demonstrated in 
Paper I that F — 1.1 works better for erg = 0.5 compared with larger values of erg. In 
any event, F = 1.1 lowers the observational estimates making the disparity with the 
observations worse. The observational estimates of WDEF are higher than ours. This 
stems from the fact that they used a different Faber-Jackson relationship for ellipticals. 
We note that the Faber et al. (1989) survey of ellipticals has even fewer bright ellipticals 
(factor ~ 1/2 for ai > 350kms~^) than we get with our Faber-Jackson relation, Mbj, = 
— 6.6364 log]^Q(cri) — 5.884. The differences are not critical since we will find, using either 
our estimates or those of WDEF, that there are too many simulated halos compared 
with the observations. WDEF study the clustering of halos with V^irc > 250kms~^ but 
they only report numbers for VcWc > 100 and 200 kms~^, so the entry for 250kms~^ is 
unfilled in Table 1. 

The triplet of numbers in Table 1 are for Vdrc > 100, 200, 250kms~^. Results arc 
shown for halos defined using EOF and DENMAX, and for DENMAX after breaking up 
the massive halos into groups and clusters of halos (see § 4.2 and 4.3). The WDEF results 
are shown before and after their special treatment of merging. Observational estimates 
are given in the first column. 
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We identify halos in our simulations with FOF linking lengths I = 0.1 and 0.2, and 
DENMAX with a 512"^ density grid. Circular velocities are defined at a comoving radius 
of 150 kpc. We find similar results using different DENMAX grids and different radius 
cuts (see Paper I). The variation in the numbers from CDM12 versus CDM13 are less 
than 10%. We see from Table 1 that the results for the DENMAX analyses are closer 
to the FOF (/=0.1) case than the FOF (/=0.2) case for the larger circular velocity cut- 
offs. WDEF reported their results before and after their prescription for restoring merged 
halos into clusters. They used FOF with an unspecified, small linking parameter and they 
defined their circular velocities using a mass within a sphere of mean density 1000 times 
the present critical density. Despite these differences, we find reasonable agreement for 
the numbers of halos with l^irc exceeding 200kms~^. However, their break-up procedure 
results in more than twice as many halos as we find for Vdrc > lOOkms"^. In any case, 
the simulations predict more than twice as many halos as there ought to be. 

In Figure 4 we show averages of the two-point correlation functions from the two 
P^M simulations for the halos with Kirc > 100, 200, and 250 kms"^ at as = 0.4. The 
mass is shown as a solid curve (also with la error bars in the bottom panel). Our 
mass correlation length agrees with WDEF: tq ~ 4 Mpc. We show the observed two- 
point correlation function {vq = 10 Mpc; logarithmic slope —1.8) as a straight solid 
line. We see a slight enhancement of the correlation length of the halos (dotted and 
dashed curves) compared with the correlation length of the mass. WDEF reported a 
similar enhancement. The two-point correlation function also has the wrong shape; in 
agreement with WDEF it is too steep on small scales and has a sag at r ~ 1.5 Mpc. 

The enhancement in ^ for the bright halos is not large enough to reconcile ug = 0.4 
with the observations. WDEF argued, and we present more arguments later, that the 
massive halos might represent clusters of galaxies that have merged. WDEF, in a compli- 
cated manner, found every halo that ever formed in the evolution of their models and then 
used a prescription for merging. (Rather than mimic their procedure in detail, we simply 
add halos to each massive system in proportion to its bound mass.) When WDEF applied 
their algorithm to their models, they found significantly more sites for galaxy formation 
compared with not breaking up the massive halos. (We list their numbers as before and 
after in the break-up column of Table 1.) They found ^ is significantly enhanced after 
break-up, and halos with Vdrc > 250kms~^ then matched the observed ^ fairly well. If 
one adds halos to massive systems, one gives extra weight to these systems which are 
more correlated than smaller systems, thereby increasing the correlation function (Kaiser 
1984). 

We show our results with the break-up of massive halos as points in the bottom panel 
of Figure 4. We show la error bars from the two simulations. This procedure introduces 
new small-scale pairs enhancing ^ on small scales. We find for one of our simulations that 
^ comes close to the observed line and that for the other simulation the enhancement at 
large r is small. We find that ro for the mass has a substantial simulation-to-simulation 
fiuctuation at r ~ 10 Mpc which is not surprising in a small box (see Figure 2). The larger 
correlation length of the halos corresponds to the simulation with the larger correlation 
length in the mass. Giving extra weight to the massive halos enhances ^ provided that 
there are significant contributions to ^ from long wavelengths. 
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Our crude break-up scheme is less ambitious than that of WDEF and it produces 
fewer hales. We include it here only to illustrate the points made by WDEF and to 
emphasize the problem associated with producing too many halos. In agreement with 
WDEF, we find that it is essential to break up the massive halos in order to approximately 
reconcile the two-point correlation function of resolved halos at erg = 0.4 with the observed 
two-point correlation function. However, if we break up the massive halos enough to 
enhance Tq to match the observations, then the resulting numbers of halos far exceed the 
observational estimates. Moreover ^(r) still has the wrong shape. These facts must be 
considered to be serious shortcomings of the model. The discrepancies must be considered 
tentative, however, because we know that a 51.2 Mpc box is too small. We use a 100 
Mpc box in subsequent sections. 

2.4- Halo Velocities 

In this subsection we present a few comments concerning the velocities of resolved 
halos. This subject is important for assessing velocity bias: the pairwise velocity dis- 
persion of the halos can much be less than that of the mass (see Carlberg, Couchman, 
& Thomas 1990). In Figure 5 we show (Tp(r) at as = 0.7 for CDM12(643,51.2,52) for 
the mass and for the halos (found with a 512^ DENMAX grid). We use the center-of- 
momentum to define the velocities of the halos. We see that Up is significantly smaller for 
the halos than for the mass. One possible source of this "velocity bias" is dynamical fric- 
tion of the halos with the surrounding medium (Carlberg, Couchman, & Thomas 1990), 
where the internal motions of particles in the halo exchange energy and momentum with 
surrounding particles or other halos. 

Much of the velocity bias, however, might actually be a statistical effect. To see 
this, in Figure 5 we also show ap for the mass with the particles from the two largest 
halos removed (dashed-dotted curve). These large halos have circular velocities (defined 
with R = 150 kpc comoving) of 981 kms~^ and 904 kms~^. We remove all of the 
DENMAX particles — even the unbound ones — which involve 7719 particles and 4603 
particles respectively. Removal of these particles is unrelated to dynamical friction. We 
see that the mass pairwise velocity dispersion is reduced significantly by this removal, and 
therefore the amount of dynamical friction required to explain the overall velocity bias is 
less than one might expect. The reason why removing the massive halos has such a large 
effect on cTp is simple (Bertschinger & Gelb 1991). The calculation of Up for the particles 
weights each pair, giving quadratically greater weight to pairs of the most massive halos. 
The pairwise velocity dispersions of these large objects are also much higher than they 
are for smaller systems. If we remove large halos we remove a large number of pairs of 
high-velocity particles. 

We now consider the important distinction between using the velocity of the maxi- 
mally bound particle from a halo (used by WDEF) and the center-of-momentum velocity. 
The short dashed curve in Figure 5 is cTp for the halos with V^circ > 192kms~^, but we use 
the velocity of the maximally bound particle; i.e. the one with the minimum potential 
computed by direct summation of particles in the halo treated in isolation. By using the 
center-of-momentum for the velocity of the halo rather than the velocity of the maximally 
bound particle, we get lower values of CTp because individual particles have a significant 
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velocity dispersion about the mean halo velocity; i.e. we are not including the internal 
halo velocity dispersion when we use the center-of-momentum. It makes sense to define 
the velocity of a halo using the center-of-momentum of the halo because observers define 
redshifts using the average velocities of the stars in a galaxy. 

3. PEAK PARTICLES 

In this section we discuss an alternative definition of "galaxies" based on particles 
initially in density peaks. This method is often used with N-body simulations lacking 
sufficient mass or force resolution to resolve evolved dark matter halos (see Kaiser 1984; 
DEFW; Bardeen et al. 1986; Park 1990; Park 1991; Katz, Quinn, k Gelb 1993). Galaxies 
are identified as particles nearest initial (linear) density maxima and their evolution is 
followed along with the other particles representing intergalactic clouds of dark matter. 
Only peaks with density exceeding some threshold are accepted. The peak threshold, for 
a given gaussian smoothing radius used to smooth the initial density field, is chosen to 
give the correct number of bright halos in the simulation volume. 

3.1. The Two-point Correlation Function: ^(r) 

Following Frenk et al. (1988), we use gaussian smoothing radii Rs of 550 kpc co- 
moving and 880 kpc comoving, with corresponding density fluctuation thresholds ly = 2.6 
and 3.0 (for Rg = 550 kpc comoving) and v = 2.0 and 2.5 (for Rg = 880 kpc comoving) 
in units of ap — the density dispersion computed from the smoothed, initial density field. 
These values of u give roughly 700 and 1600 galaxies, respectively, in a (100 Mpc)^ vol- 
ume (our simulation CDM16). Using the parameters in the Schechter luminosity function 
(see Paper I where we used parameters from Efstathiou, Ellis, & Peterson 1988) these 
correspond to circular velocity cut-offs of approximately 250 km s"-*^ and 200 km s"-*^ re- 
spectively. For a given gaussian smoothing radius smaller values of u correspond the 
smaller circular velocity cutoffs. We choose two values of Rg to test the sensitivity of our 
results to this parameter. 

We next determine ^ computed using only the peak particles in the CDM16 simu- 
lation. In other words, we compute the two-point correlation function using the present 
positions of the particles which are tagged as galaxies. The results are shown in Figure 6. 
The peak particles indicate that ag = 0.5 is possibly suitable as the present epoch — the 
value of To is roughly 10 Mpc and the logarithmic slope is very nearly —1.8 from about 
1 Mpc to 20 Mpc. Note that this success is exactly what led DEFW to champion biased 
CDM. The correlation lengths for erg = 0.7 and erg = 1.0 are also roughly 10 Mpc, but the 
slope steepens at roughly r ^ 3 Mpc for ug = 0.7 and at roughly r ;S 4 Mpc for ug = 1.0. 
Even for ug = 0.5, the slope steepens for roughly r ^ 1.5 Mpc. (We also computed ^ 
at (Tg = 0.4 and the results are nearly identical to erg = 0.5 except that the steepening 
of the slope occurs at r ^ 1.25 Mpc rather than r ~ 1.5 Mpc.) The enhancement occurs 
because peak particles are more likely to be found in massive halos where the chance of 
a peak being above the threshold u is higher. 

To see how peaks are associated with massive halos, we show in Figure 7 the bound 
particles from a massive halo (2.1 x 10"^^ M©) at ug = 0.5 from CDM16 (upper left 
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panel). We use a 512^ DENMAX grid which apparently has not completely resolved all 
substructure. We noted this problem in Paper I where we concluded that increased force 
resolution reveals substructure and increased DENMAX grids are required to bring out 
this substructure. However, in many cases there is no obvious substructure in the images 
of the massive halos. In the upper right panel we show peak particles that end up as 
bound members of the massive halo shown in the upper left panel. We see that there are 
many peak particles in this massive halo. (We discuss the other panels later.) 

The large number of peak particles per massive halo is typical. Conservation of 
numbers then implies that peak particles must undersample less massive halos outside 
clusters. To see this quantitatively, we consider the resolved halos found by DENMAX 
with Fcirc > 250kms~^ defined at i? = 200 kpc comoving at erg = 0.5. There are 737 
halos. Of these, we count the number of halos that do not contain any peak particles as 
bound members. For Rg = 550 kpc comoving, using u = 3.0 which yields 639 peaks, 425 
resolved halos have no peak particles as bound members. (Using u — 2.9 yields 826 peaks, 
and then 362 resolved halos have no peak particles as bound members.) For Rs = 880 
kpc comoving, using u = 2.5 which yields 740 peaks, 355 resolved halos have no peak 
particles as bound members. We conclude that about half of the massive halos in the 
evolved, nonlinear density field contain no peak particles. This is a major failing of the 
peak particles as galaxy tracers (cf. Katz et al. 1993) and it calls into question N-body 
simulations that rely on peak particles in lieu of dense halos. 

3.2. Discussion 

Should we conclude from Figure 6 that the correlation function slope in the CDM 
model is too steep on small scales? Park (1991) presented similar studies of ^(r) using 
peak particles, but he did not show the steepening of ^ within 1 Mpc. His force resolution 
was of order ~ 1 Mpc (for a 256^ grid PM simulation in a 153.6 Mpc box). However, his 
models show, in agreement with our results, the steepening of the slope at larger scales 
for CTg = 1.0. We will have to compute ^(r) using actual resolved halos, with and without 
the break-up of merged massive halos, to decide whether the steepening of ^ (r) is real or 
an artifact of peak particles. 

In any case, we conclude that the method of peak particles can give misleading 
results. In § 2 we found that we could only get a significant enhancement in ^ at erg = 0.4 
if we broke up massive halos, but this produced far too many halos. We found that we can 
get the required enhancement in ^ at ag = 0.5 (the results were similar at erg = 0.4) using 
the correct number of peak particles. However, half of the actual, massive, nonlinear halos 
did not contain any peak particles. Because there is nothing unusual about those halos 
that do not contain peak particles (compared with halos of the same circular velocity that 
do), we cannot argue that peak particles are to be preferred over direct identification of 
dense resolved halos. Peak particles oversample the clusters and undersample the field. 
These effects enhance the two-point correlation function with fewer halos. For these 
reasons, we must study the CDM model using resolved dense halos, rather than peak 
particles, to trace galaxies. 

4. TWO-POINT CORRELATIONS OF HALOS 
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4.1- Introduction 

The preceding section has motivated a study of the two-point correlation function 
using resolved halos, which we undertake in this section. We use the high-resolution 
simulation CDM16(144"^,100,85) in an attempt to constrain erg based on the slope and 
amplitude of ^(r). We will find that the results depend on how merged halos are treated, 
so we will devote some discussion to this issue. 

We show ^ computed from resolved halos at ag = 0.5, 0.7 and 1.0 from CDM16 in 
Figure 8. We use a 512^ DENMAX grid, with bound particles only (see Paper I), for 
the remainder of this paper. We see that the correlation length, tq, falls short of the 
observed value at low values of cxg. Also, there is increased merging at later epochs (see 
Paper I); this explains why the halos are antibiased, i.e. are less clustered than the mass, 
on small scales. The antibiasing is stronger at later epochs and for smaller halos. This 
is because merging increases with increasing erg and the smallest halos merge into larger 
systems. Unfortunately for us, observers do not directly measure the clustering of the 
mass; they measure the clustering of the galaxies. However, we see that unless galaxies 
are clustered more strongly than the halos, we will not be able to match the observed 
two-point correlation function. 

Carlberg & Couchman (1989) performed a simulation with 1.2 x IO^^Mq particles 
in a 80 Mpc box. (CDM16 has particles with 2.3 x lO^^M© in a 100 Mpc box.) At 
CTg = 0.54, using FOF to identify dark halos, they found, as we do, that the dark halos 
are antibiased with respect to the mass on small scales and that they trace the mass on 
larger scales (see their figure 8b). Couchman & Carlberg (1992) studied more evolved 
models and they also found the same level of antibiasing of the dark halos with respect 
to the mass on small scales. 

If we take these results for ^ at face value, then the = 1 CDM model has seri- 
ous shortcomings: the correlation length is too small for ag ~ 0.7; and the correlation 
amplitude is too small and turns over on small scales, particularly for ug ^ 0.7. Rather 
than abandon the model, however, we explore the possibility that restoring merged halos 
might sufficiently increase ^. This step is reasonable, because the most massive halos 
clearly ought to contain several galaxies. However, we will find that we create as many 
problems in the process as we solve! 

We break up the massive halos using two techniques. The first method (fall-in 
method) involves finding the maximally bound particle in each halo at an earlier epoch. 
We identify the maximally bound particles that have fallen into massive halos represent- 
ing clusters at the present epoch. The second method (mass-to-light method) involves 
assuming a mass-to-light ratio for the massive halos representing clusters, and then, using 
the observed Schechter luminosity function, assigning the appropriate number of mem- 
bers to the clusters. The fall-in method directly shows that smaller halos merge into 
larger halos, as we showed in Paper I. Perhaps dissipative effects (Katz et al. 1992; Katz 
& White 1993; Evrard ct al. 1994) or harder forces (Dubinski & Carlberg 1991) might 
help these systems survive the merging process. However, because clusters exist but are 
over-merged in our simulations, it is reasonable to unmerge the most massive halos. 

4-2. Break-up of Halos: Fall-in Method 
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In the fall-in method we find all of the bound 512^ DENMAX objects at an early 
epoch, which we call the tagging era, and we find the maximally bound particle from each 
of these halos. We then grab the present day positions and velocities of these particles and 
we add each one to our list of present day halos only if it is a bound member of a massive 
halo with present day circular velocity (defined within a radius 200 kpc) > 350 kms"-*^. 
Thus, wc break up the bound mass of large halos into several different objects that were 
distinct entities at the tagging era. Each such object is assigned the circular velocity it 
had at the tagging era. We also retain the massive merged halo, unless the sum of masses 
of the added halos exceeds its mass. 

There are two arbitrary parameters in the method: the circular velocity beyond 
which we break up the massive halos and the epoch which we choose as the tagging 
era. For the former we choose Vcirc = 350 km s"-*^. For larger values we found excessive 
numbers of halos in Paper I. We are admittedly forcing improved agreement with the 
observations. For the tagging era we try erg = 0.2 and erg = 0.3. There is no ideal, single 
epoch since galaxy formation is a continual process. WDEF eliminated this ambiguity 
by finding every halo that ever formed and then by putting merging in by hand. We do 
not attempt to reproduce their procedure. 

We show an example of the fall- in method at ag = 0.5 in the lower left panel of 
Figure 7. This example uses erg = 0.2 for the tagging era. Ten objects with V^irc > 
250 km s-^ fall into the massive halo shown in the upper left panel. When we use erg = 0.3 
for the tagging era there are four halos with V^circ > 250kms~^ that fall into this massive 
halo. Larger values of ug contain larger objects, but many of these objects have already 
undergone merging. 

We now examine the effect of the fall-in method on ^ at erg = 0.5. We can see in 
Figure 7 that we introduce more pairs on small scales (lower left panel) and that we give 
extra weight to the massive halo (upper left panel). We show ^ after applying the fall-in 
method to CDM16 in the top panel of Figure 9. (The M/£, method is discussed in the 
next section.) Compare these results with the results without break-up in the top panel 
of Figure 8. We see better agreement of the slope with the observed slope except on scales 
< 1 Mpc. We also see that we increase the correlation length closer to the observed value. 

The numbers of halos with no special treatment of the massive halos at ug = 0.5 
are 1340 for Kirc > 200 km s"^ and 737 for Kirc > 250 km 8"^ The numbers from the 
fall-in method at erg = 0.5 are 1934 (1706) for Kirc > 200 km s"^ and 1022 (940) for 

Vcirc > 250kms~^ if we use erg = 0.2 (erg = 0.3) as the tagging era. There is a ~ 30 — 40% 
increase in the numbers of halos using the fall-in method, showing that there has been a 
significant amount of merging. 

The fall-in method and the peak particles method produce nearly the same shape 
of ^ (see Figure 6). However, the steepening of the slope occurs at r ~ 0.8 Mpc using 
the fall- in method, compared with r ~ 1.5 Mpc using the peak particles method. Both 
produce the wrong shape. 

Notice that the feature at r ~ 1.5 Mpc in the bottom panel of Figure 4 is more 
prominent than it is in Figure 9. This is because the halos are more extended using 
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a 512'^ DENMAX analysis in a 100 Mpc box compared with using a 512^ DENMAX 
analysis in a 51.2 Mpc box, as we discussed in Paper I. The effectively coarser DENMAX 
allows more peripheral particles. These peripheral particles make the halos bigger and 
introduce more pairs beyond ~ 1.5 Mpc. 

The most important difference between the fall-in method and the peak particles 
method is that the former requires far more halos to get the same level of enhancement 
as the peak particle method. This is to be expected because, as we demonstrated in § 
3, the method of peak particles oversamples the clusters and misses many field galaxies. 
If galaxies cluster like dense dark matter halos, then it would seem that the amount of 
clustering bias (i.e., the ratio of ^ for galaxies to that for the mass) must be less than 
predicted based on peak particles. 



As an alternative to the fall-in method for breaking up massive halos we consider 
an ad hoc method designed to constrain the mass-to-light ratio Mj L of clusters. We 
associate these galaxy clusters with the massive merged halos and assign each such halo 
the number of galaxies expected on average given a universal luminosity function. This 
method sacrifices all predictive power for cluster M/£'s, but we do not consider this a 
grave loss because we doubt that any reasonable attempt can be made to estimate the 
luminosities of galaxies in a cluster using a purely dissipationless simulation that follows 
only the dark matter. With the Mj C method, we assume only that the most massive 
halos should be associated with galaxy clusters and that a specified Mj L applies for all 
such clusters. This simple-minded prescription offers, at least, a useful foil for the fall-in 
method. Moreover, it allows us to vary the richness of clusters by varying a single number, 
Mj C On the other hand. Ashman, Salucci, & Persic (1993) argue that observations of 
disk galaxies imply a variable Mj C which could reduce the excess richness and numbers 
of galaxies in clusters implied by hierarchical models assuming constant Mj C 

Our procedure is straightforward. We examine all halos with V^xyc (defined at = 
200 kpc comoving) exceeding 350kms~^. For each such halo we divide its total bound 
mass by a specified Mj L (with L measured in the blue) to get the total blue luminosity in 
the cluster. Ramella, Geller, & Huchra (1989) find MjL ~ 180/i (in units of M©//:©) for 
groups in the CfA2 survey. Some clusters are estimated to have values exceeding 500/i, 
but there is still controversy among workers in the field. Trimble (1987) gives a review. 

We obtain a distribution of galaxies using the Schechter luminosity function $(>C) 



with parameters = 1.56 x 10~^/i^ Mpc"^, = -19.68 - 2.51ogio/i"^, and a = 

— 1.07 (Efstathiou, Ellis, & Peterson 1988). The total luminosity in a volume V is 



4-3. Break-up of Halos: Mass-to-light Method 



oo 




(4.1) 







The total number of galaxies in a volume V with a luminosity exceeding £ is: 



oo 




(4.2) 
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Combining eq. (4.1) and (4.2) and defining x = C/C^, we get the total number of halos 
exceeding a luminosity £ in a cluster with total light >Ctotai: 

JV(>£,A„.aO = ^%^. (4.3) 

Colless (1989) and Schechter (1976) inform us that the same luminosity function works 
for rich clusters and for field galaxies within the uncertainty of the data. 

We now put the steps together. We take the bound mass of one massive halo (those 
with Fcirc > 350kms~^) and we divide it by a specified universal value of M/C This 
gives us the total luminosity emitted by the cluster: £totai- We then add A^(> £, £totai) 
halos with luminosity exceeding C to the big halo using eq. (4.3). We relate C to circular 
velocity using the TuUy-Fisher (Pierce & TuUy 1988) and Faber-Jackson relations (using 
our fit from Faber et al. 1989; see paper I for details). We assume 70% spirals and 
30% ellipticals; for the latter, the Fcirc is corrected to ai using eq. (2.2) with F = 1 (no 
significant difference occurs if we use F = 1.1). The value of C corresponding to ^circ is 
chosen so that the number in eq. (4.3) exceeding C is the same as the number exceeding 

When we add in halos using this mass-to-light method we need to choose positions 
and velocities. We do this by randomly sampling the massive halos we are breaking up. 
In other words, if the massive halo contains Nh particles, we generate a uniform random 
number from 1 to and we use the present day position and velocity of particle A^^- 
We repeat this procedure for each added halo. The break-up of the massive halo in the 
upper left panel of Figure 7 is shown for the mass-to-light method with random position 
sampling for M/C = 125 in the lower right panel of Figure 7. In Figure 9 wc show ^ at 
CTg = 0.5 for M/C = 125 where we did this random sampling (long dashed curve), and 
where we put all added halos on top of each other at the locations of the original massive 
halos (dot-dashed curve). The results agree at larger scales, but it is essential to use the 
random sampling method to see the efifects from close pairs. We use random sampling 
for the remainder of this paper. 

In Figure 9 we see that the mass-to-light method produces results similar to the fall-in 
method and the peak particles method. We notice, however, that the slope on small scales 
is steeper for the mass-to-light method than for the fall-in method but comparable to the 
peak particles method. This is because there are more galaxies added with M/C = 125 
than with the fall-in method. Wc quantify these numbers later. By varying M/C we can 
test whether the spatial and velocity statistics as well as numbers of halos are acceptable 
for a given model. These quantities are not guaranteed to all work out satisfactorily even 
with the freedom we allow ourselves in how massive halos are broken up. We will see, 
on the contrary, that small-scale galaxy clustering and velocities, combined with galaxy 
abundances and group multiplicities, present serious difficulties for the CDM model. 

4.4- Constraining as Using ^{r) 

In this subsection we investigate the two-point correlation function of simulated 
galaxies from our 100 Mpc high-resolution CDM simulation CDM16(144^, 100,85). Galax- 
ies are identified with DENMAX halos except for the massive halos, which are split into 
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several galaxies using either the fall-in method or the mass-to-light method, with Mj L 
a parameter that we can vary. The purpose is to determine whether there exists a nor- 
malization (T8 such that ^ for the simulated galaxies matches the observations. We vary 
the break-up procedure to determine the sensitivity of our conclusions to this uncertain 
step. 

Figure 10 shows results at three epochs erg = 0.5, 0.7, and 1.0. The shape of ^ fails to 
match the observed ^ in all cases and this must be considered to be a serious shortcoming 
of the models. At a% = 0.5 the enhancement in ^ is nearly sufficient for M/C = 50 but 
the slope is too steep on small scales. Gott & Turner (1979) showed that the logarithmic 
slope —1.8 is valid down to at least scales of ~ 10 kpc with no indications of any features 
on small scales. The M/C = 250 case at ug = 0.5 is not too steep on small scales, but 
the correlation length is only ^ 7 Mpc. For M/C = 500 the correlation length is ~ 6 
Mpc and ^ falls between the no break-up case and the M/C = 250 case at small scales 
at (Tg — 0.5. The no break-up case at (Tg = 1 is almost acceptable, but the significant 
turnover on small scales does not match the observed slope and the massive halos do 
not look anything like observed clusters, i.e. they are single objects rather than tens of 
objects. At no epoch, with no treatment of halo break-up, do the simulations match the 
observations. 

We now examine the numbers of halos and the properties of the clusters since these 
are central to further conclusions regarding ^. In Table 2 we list the numbers of halos 
with Vcirc ^ 250kms~^. The numbers are for the (100 Mpc)^ volume. The numbers 
without break-up are in the default column; we also show numbers for M/C = 50, 125, 
250, and 500. 

The observed number is less than 621 (563) for VcWc > 250kms~^, assuming ai = 
Vcirc/iV^/F) for ellipticals with F = 1 {F = 1.1). (Again, this is an overestimate 
because of the assumed Faber- Jackson relationship; sec Paper I.) Even before breaking 
up the merged halos the numbers are too large; breaking up the halos leads to an even 
greater disagreement with observations. The numbers for M/C — 50 and 125 are factors 
^3 — 10 too high! Therefore, unless the mean M/C ^ 1000/j. for typical groups, we can 
safely rule out ug ^ 0.7 just from the numbers shown in Table 2. 

The numbers for M/C = 250 at erg = 0.5 are also high, but may be consistent with 
observations within various uncertainties. If we choose M/ C ^ 250 we partially solve the 
high galaxy abundance problem and the correlation function steepness problem, but we 
do not raise the correlation length to the observed value. Remember, WDEF found a 
factor of ~ 3 too many halos to yield the required enhancement in ^ at erg = 0.4 (although 
they used a 50 Mpc box), and we see in Figure 10 that the correlation length falls short of 
the observed value at erg =0.5 for M/C ~ 250. At later epochs we can solve the steepness 
problem using the catalogs without breaking up the clusters, but then our simulations 
do not have rich clusters like the real universe. 

Finally, the numbers of galaxies from the fall-in method at erg = 0.5 are comparable 
to the M/C — 250 numbers in Table 2 at erg = 0.5. This explains why ^ at erg = 0.5 
looks markedly similar for the fall-in method and for the mass-to-light method with 
M/C = 250. This lends some support to our use of the mass-to-light method. If gaseous 
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dissipation is able to preserve galaxies in clusters, even when the dark matter halos merge 
(White & Rees 1978; Katz & White 1993; Evrard et al. 1994), the CDM model might 
produce clusters of galaxies with Mj L 250, in not too violent disagreement with the 
observations. 

5. CLUSTERS AND THEIR RICHNESS 

It is not enough for CDM or any other theory to predict the correct two-point 
correlation function, pairwise velocity dispersions, and galaxy abundances. A successful 
theory must also predict the correct richnesses, abundances, and mass-to- light ratios of 
galaxy groups and clusters. As we have noted, this test is difficult to make using a 
purely gravitational N-body simulation without dissipation because of the overmerging 
problem. However, we showed in the preceding section that a plausible scheme for undoing 
the overmerging is based on assigning galaxies to massive merged halos in proportion to 
the mass. Is it possible to do this with a reasonable value of Mj C so that the correct 
group multiplicity function (richness distribution) is obtained? Are there then too many 
clusters? 

Let us recall first that the mean Mj L for = 1 is ~ 750 for h — 0.5. On the 
other hand, most dynamical measurements on cluster scales yield much values smaller 
by a factor of three or more (see, e.g., Peebles 1986). While velocity bias in clusters 
might reduce the apparent Mj L to acceptable values for = 1, there exist more direct 
mass measurements from X-ray emission for some clusters (e.g., Hughes 1989). We will 
therefore examine as well how much mass is contained in our massive halos. 

We compute the fraction of the total mass in our (100 Mpc)*^ volume contained in 
massive halos at erg = 0.5, 0.7, and 1.0 using CDM16. We accumulate the total bound 
mass in all halos with V^irc (defined at 200 kpc comoving) exceeding 350kms~^; these 
are the objects that we have been breaking up in the previous section. We find the 
percentage of mass contained in these objects to be 19.2% (267 objects) at o% — 0.5; 
29.9% at 0-8 = 0.7 (363 objects); and 39.9% (420 objects) at ag = 1.0. These numbers 
are excessive when one recalls that only a few percent of galaxies are in rich clusters; see 
Bahcall (1979) for a review. 

Since these fractions are so high, we compute a few more interesting numbers. We 
compute the fraction of the mass contained in objects at = 1.0 using larger circular 
velocity cut-offs. For objects with V^irc > 400kms~^ the mass fraction is 36.9% (301 
objects). For objects with V^irc > 500kms~^ the mass fraction is 31.5% (170 objects). 
Therefore, the amount of mass contained in very massive objects is enormously high. In 
Paper I we learned that the cumulative mass fraction converged with increasing mass 
resolution if we imposed a distance cut. Therefore, we compute the mass fraction of 
objects above a given circular velocity cut-off defined at 500 kpc comoving, and we 
accumulate only the bound mass within 500 kpc comoving. At erg = 1 for V^irc > 
350kms~^ the mass fraction is 20.2% (344 objects). At cxg = 1 for Vaxc > 500kms~^ the 
mass fraction is 15.0% (156 objects). The numbers of objects are slightly less for cut-off 
radii of 500 kpc comoving versus 200 kpc comoving because many of the circular velocity 
profiles with V^irc ~ 500kms~^ are actually falling slightly at these scales. However, for 
larger circular velocities the profiles are still rising beyond 200 kpc comoving. 
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There is some uncertainty in defining the bound mass of our massive halos. However, 

even using a conservative estimate wc find that at least 15% of the mass is contained 
in very massive halos (with V^circ > 500kms~^) at erg = 1. On the other hand, the 
percentage is not 100% so we can consider mass-to-light ratios < 750 for our massive 
objects (remember that f2 = 1 demands Mj L = 750 on average) if Mj L > 750 for less 
massive objects. Unfortunately for Q = 1 CDM, this goes against observations (Trimble 
1987). One cannot argue that the missing mass is far outside galaxies in the CDM model 
(at least with ug ^ 0.5) because more than half the mass is within 500 kpc from the 
center of a halo (cf. Paper I). 

Next, we consider the richness of our hand-made clusters and we impose further 
constraints on the mass-to-light ratios; the reader is reminded that the numbers of halos 
in our volume also impose constraints (see Table 2). 

Ramella, Geller, & Huchra (1989; hereafter RGH) studied groups of galaxies from 
the -B(O) < 15.5 CfA2 redshift survey. Redshift space projection effects may bias some 
group properties relative to groups selected using distance information as we do in our 
simulations (Nolthenius & White 1987). Nevertheless, we feel it is useful to present 
a preliminary analysis of the group multiplicity function in the CDM model. For our 
discussion in this section we convert all relevant quantities to Zwicky magnitudes using 
B{0) ^ Bt + 0.29 (Efstathiou, Ellis, & Peterson 1988). We replicate our (100 Mpc)^ 
volume using periodic boundary conditions into a (250 Mpc)^ volume. We then select 
a wedge corresponding to the CfA2 sky coverage: right ascension range 8^^ < a < 17^^ 
and declination range 26.5° < d < 38.5°. We refer to this as the 12° slice. We assume 
Hq = 50kms~^ Mpc~"^ and we impose a distance cut of R < 240 Mpc in our analysis. We 
use actual galaxy positions rather than redshifts and we impose an apparent magnitude 
limit of -B(O) < 15.5. We assume a Tully-Fisher relationship, see Paper I, converted to 
Mb{o)i where we determine the circular velocities of the halos from CDM16 at 200 kpc 
comoving. 

We use DENMAX to identify aU halos with Kirc > V^l^ = 50kms-^; then we 
use EOF to identify groups of halos in our wedge after breaking up the massive halos 
{Voire > 350 kms~^) using the mass-to- light method. We determine a EOF linking length, 
I in Mpc, corresponding to a given galaxy overdensity Sp/p given by = 2/{nSp/p) (see, 
for example, Frenk et al. 1988) where n is the number density of halos with circular 
velocity exceeding V^^™ from our original (100 Mpc)^ volume. We use EOF to identify 
groups of halos after breaking up the massive halos, but prior to imposing an apparent 
magnitude limit. Typical values of Z, for (5p/p = 80, range from 0.8 Mpc to 1 Mpc for the 
various assumed values of M/C and erg. 

We only identify groups with three or more members to be consistent with RGH. 
RGH chose a linking distance using a galaxy number density estimated from the observed 
Schcchter luminosity function. However, they varied their linking length with redshift to 
account for the sparse sampling of galaxies at large redshift. We avoid this difficulty by 
applying EOF with a fixed linking length prior to applying an apparent magnitude limit. 
We then apply the apparent magnitude limit to the resulting group catalog in a manner 
described below. 
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For field halos, i.e. those that are not in groups with 3 or more members, we 
simply compute Mb(o) using the TuUy-Fishcr relationship, and we remove those with 
-B(O) > 15.5. For the halos in groups we apply the following procedure. If the group 
member is not created from the break-up of a massive halo, then we eliminate it if 
B{0) > 15.5. For group members that are created from the break-up of a massive halo, 
we remove all of them and replace them by the number of halos determined from eq. (4.3) 
for an assumed, universal M/C The lower luminosity limit in eq. (4.3) is computed from 
15.5 — Mb{o) = blogigd + 25.0, where d is the distance to the group centroid in Mpc. 
Note that here we do not need to relate luminosity to V^irc in eq. (4.3). However, to be 
consistent with our use of V^^^^, we never allow C to fall below >Cmin determined from 
V^l^ using the TuUy-Fisher relationship. 

The basic parameters in the group finding algorithm are the galaxy overdensity Sp/p 
used to determine the linking parameter, the faint cut V^^™, the mass-to-light ratio M/C 
used to break up the massive halos, and the circular velocity cut-off above which we break 
up massive halos. We discuss these four parameters here. 

1) We report results using 6p/p = 80, the middle value considered by RGH, since we 
see the same levels of variation with 6p/p as reported by RGH and our conclusions do 
not depend critically on its value. 

2) We report results for — 50kms~^. Our results do not depend sensitively 
on V^l^ because the low mass galaxies quickly fall out of sight. For example, in a case 
where we identify 1555 field galaxies in our 12° slice with an apparent magnitude limit, 
only 233 have Vdrc < 125kms~^ and only 62 have Vdj-c < 75kms~-^. This is encouraging 
since we found in Paper I that we had factors ~ 2 — 3 too many halos compared with 
the observations for V^irc ^ 125kms~^. In a magnitude-limited survey we would not be 
swamped by low mass halos. 

3) We report results using M/i2=125, 250, and 500. From a list of 36 groups, RGH 
found a median M/C of 178/i = 89 for h = 0.5. We choose large values of M/C because, 
as we will see, even M/C=125 produces groups that are too rich. 

4) There is some arbitrariness to the value of Vdrc above which we break up the 
massive halos. We report results using V^irc = 350kms~-^. If we raise this value we get 

too many isolated massive halos (see Paper I) which are too big to represent individual 
galaxies. On the other hand, the numbers of halos added quickly approaches zero below 
Vciic = 350kms~^ for the M/C studied here. 

The results from our simulations are shown in Table 3. We report numbers from 
RGH for the full 12° slice, but we impose a redshift cut of 12000 km s~^. RGH only 
studied groups with centroids < 12000 km s~^. We report numbers from the simulations 
for the full 12° slice for R < 240 Mpc. The table shows the number of groups, iVgroups, 
identified with 3 or more members, with 10 or more members, and with 20 or more 
members. We also show the number of galaxies, A^gaiaxies, in the field, i.e. those that are 
not in groups with 3 or more members. We estimate the number of CfA2 field galaxies 
within 12000 kms~^ as follows. The CfA2 catalog has 1766 galaxies and we estimate 
from figure 1 in RGH that ^ 100 galaxies are beyond 12000 kms-^ RGH found 778 
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galaxies in groups with three or more members and only a handful of these galaxies are 
beyond 12000 kms~^. Therefore, the number of field galaxies within 12000 kms~^ in 
the CfA2 catalog is approximately 1766 — 778 — 100 ~ 900 galaxies. The last entry in 
the table, is a richness statistic defined later. The reader is cautioned that RGH 

estimate that ^ 30% of the groups with 3 or 4 members might be an artifact of projection 
effects. 

We can draw several important conclusions from the results shown in Table 3. If we 
do not break up the massive halos, then we do not have enough groups and there are no 
groups with 10 or more members. Therefore we need to break up our massive halos if our 
simulated universe is to contain groups comparable to the observed numbers! In all cases 
we have too many field galaxies. We demonstrated earlier that these are not dominated by 
faint galaxies. However, in Paper I we found that we had the correct number of halos with 
circular velocities between 150kms~^ and 350kms~^. The reason for this discrepancy is 
that here we apply only the TuUy-Fishcr relationship to the halos (i.e., we are treating 
all halos as spirals) rather than a combination of the TuUy-Fishcr relationship and the 
Faber-Jackson relationship as we did in Paper I. Applying the TuUy-Fisher relationship 
to elliptical galaxies, which tend to be the most massive halos, makes the halos appear 
brighter than they really are. On the other hand, most of our group members result from 
the break-up of massive halos where we do not need to assume a relationship between 
circular velocity and luminosity. Because of this problem, we should give more emphasis 
to the richness of our groups than to the apparent excess of field galaxies. (We use the 
TuUy-Fisher relation only for field galaxies; massive halos are broken into galaxies based 
on an assumed mass-to-light ratio.) 

We can constrain Mj L ^ 125 = 250/i based on the number of groups with 3 or more 
members and the total number of galaxies in all groups with 3 or more members. For 
CTg = 1, M/C must be ~ 250 = bOOh. In most cases, however, we still have too many rich 
groups with 50 or more members. We should note that because the observed number of 
groups with three or four members may be contaminated by projection effects, the total 
numbers of objects in groups could be smaller by ^ 30% (see RGH). This would lower 
the observed numbers in groups and, by definition, raises the observed numbers in the 
field, although it does not solve the problem of too many rich groups. 

To further quantify the richness of our groups, we compare the cumulative number 
of galaxies in groups with the estimates from RGH for the CfA2 survey. The cumulative 
number of galaxies in groups is defined by RGH as: 

— -^m e mb e rs 

-^galaxies 

N=3 

where A''gaiaxies(^ A^members) IS the total number of galaxies contained in all groups with 
three to A^members mcmbcrs and Ng{N) is the number of groups containing A^ members. 
The results are shown in Figure 11 using M/£ = 125, 250 and 500, and the no break 
up cases, at ug = 0.5, 0.7, and 1.0. Figure 11 is computed for a 6° slice (we divide the 
numbers from our 12° slice by two) to compare with RGH using a 6° (their figure 2). 
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We clearly see the dramatic shortcoming of the no break up cases at all epochs. 
When the massive halos are broken up we find groups that are richer than observed by 
RGH; the rise in the predicted cumulative galaxy number is also generally slower than 
the results for the CfA2 survey indicating that our group members are concentrated in 
relatively larger groups. A useful statistic is shown in Table 3. This is the value of 
-^members whcre the cumulatlvc number of galaxies in groups reaches 1/2 its maximum 
value. The value of indicates that we need M/ C Zj 250. We can rule out M/ C = 125. 
The remaining question is whether or not nature can hide a lot of mass; this will be an 
important consideration when we study velocities in the next section. 

6. HALO PAIR VELOCITY DISPERSIONS AND CLUSTERS 

We now consider constraints on erg from CDM16(144^,100,85) based on pairwise 
velocity dispersions of the resolved halos. We address the following questions: What is 
CTp for the halos without the break-up of massive halos? How do we assign velocities 
to halos added to the massive halos, and what is the effect on Up? Is there a linear 
normalization of the = 1 CDM model, (Tg, when the pairwise velocity dispersions agree 
with the observations? 
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6.1. Constraining Using a ^{r) of Simulated Halos 

The pairwise velocity dispersions of the halos from CDM16 without breaking up the 
massive halos are shown in Figure 12 at ug = 0.5, 0.7, and 1.0. We define the circular 
velocities at 200 kpc comoving, and we show results for VcWc > 100, 150, 192 and 250 
kms^"*^. 

The open symbols are the observed estimates from the Davis & Peebles (1983, here- 
after DP) analysis of the CfA B{0) < 14.5 redshift survey. The different symbols are for 
different modeling parameters. The best estimates are open circles with la error bars 
(shown as vertical lines). The squares are for a different set of modeling parameters, and 
the triangles are results with three clusters removed. The details are not important for 
our purposes; the scatter is small compared with the erg dependence of fXp. The results at 
r ~ 10 Mpc are the least accurate because of distortions from peculiar motions. CDM16 
has a Plummer softening of 65 kpc comoving which affects small scales. For these reasons, 
we will focus our comparisons at r ~ 1, 2.2, and 4.6 Mpc. 

DP removed all galaxies with Mb{o) > —18.5 + 51ogiQ h = —20. If we convert this 
to the Bt system and use the TuUy-Fisher relationship (see Paper I), this corresponds to 
removing all halos with Vdrc ^ 175kms~-^. We study all halos with Vdrc > 150 km s"-*^ 
and Vcirc ^ 250kms~^. The former is important since the pairwise velocity dispersions 
increase with increasing circular velocity cut-off and simulated dispersions are higher than 
the observed estimates at ^ 0.7. If we are to rule out any values of erg, it is better to 
be conservative. 

Based on Figure 12, observational data constrains erg 0.7. The case ug = 0.5 is an 
excellent match to the observed data. The results are in reasonable agreement with the 
observed data at ag = 0.7 for Vdrc ^ 100 and marginally for V^irc > 150kms~^. The case 
(Tg = 1.0 is ruled out; the pairwise velocity dispersions are too high by factors ~ 1.5 for 
r ^ 1 Mpc. Note that this is true even though there is a velocity bias of about a factor 
of two! 

There are two important issues we need to consider. We notice that the pairwise 
velocity dispersions are significantly lower for the resolved halos than for the mass. This 
velocity bias was discussed in § 2.4. We also need to compare our results with Couchman 
& Carlberg (1992, hereafter CC) who investigated erg ^ 1.0 with a 2 million particle 
P^M simulation in a 200 Mpc box. CC used a different definition to normalize the linear 
CDM power spectrum; their bi — 0.8 corresponds to ug ~ 1.0. CC assumed Q — 1, 
Ho = 50kms~^ Mpc~"^, and their particle mass is 2.65 x 10^^ M© compared with our 
particle mass of 2.3 x 10^° Mq for CDM16. 

CC found a pairwise velocity dispersion for halos with M ^ 2.1 x lO^^M© at 1 Mpc 
of ~ 490kms~^ in agreement with our results for the lower circular velocity cut-offs. CC 
found a pairwise velocity dispersion for the mass of ~ 2300kms-Vv/3 ^ 1325 km s'^ at 
1 Mpc; this is again in agreement with our results. CC did not report pairwise velocity 
dispersions on larger scales where we find the disparity with the observations to be large. 
CC also found that their halos have smaller two-point correlations than the mass; this is 
in agreement with our results presented in § 4. 
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We argued in § 4 that we need to break up massive halos into clusters to remove the 
turnover of ^ on small scales and to enhance the correlation length, and in § 5 because 
clusters really exist in our universe. We also argued in § 2.4 that using the center- 
of-momentum of resolved halos significantly reduced the pairwise velocity dispersions 
compared with the mass since a significant number of high velocity particles are contained 
in a few massive halos. For these reasons it is important to consider the effect on Up of 
breaking up the massive halos before further conclusions can be drawn. 

We use the mass-to-light method to break up the halos with Vdrc > 350 km s""*^, 
randomly sampling the positions and the velocities of the massive halos to assign positions 
and velocities to the added halos. The results are shown in Figure 13. We see immediately 
that the pairwise velocity dispersion for the halos now traces that for the mass. Wc 
have introduced a significant number of pairs with high velocity dispersions; the added 
cluster members sample massive halos which have high velocity dispersions. These results 
indicate that the pairwise velocity dispersions are too high at ag = 0.5, 0.7, and 1.0 if we 
break up the massive halos. 

CC did not report the high pairwise velocity dispersions associated with clusters. 
They found that merging decreases the numbers of halos in high dispersion regions, and 
they referenced Bertschinger & Gelb (1991) where we first discussed why this effect can 
significantly reduce pairwise velocity dispersions. However, CC did include a prescription 
for preserving merged systems as distinct halos found by FOF, but they commented that 
only ~ 20% of their "galaxy precursors" survive as distinct "galaxies" . A group analysis 
of the CC data, as we have done in § 5, is needed to estimate their group multiplicity 
function. Because our default catalogs (no break-up) reveal Up in agreement with CC at 1 
Mpc, we suspect that they would see higher ap if they had the requisite group multiplicity 
function. 

Before we can rule out any values of erg we must examine lower circular velocity 
cut-offs. We must also consider the possibility that the velocity dispersions of galaxies 
in clusters can be less than the velocity dispersions of the dark matter. Finally, we must 
consider M/C — 500 at as — 0.5 which compares favorably with the observed properties 
of groups of galaxies. These tests are the focus of the next subsection. 

6.2. Velocities of Added Cluster Members 

In the previous section we randomly sampled the velocities of the particles in the 
massive halos to assign velocities to the added cluster members. An alternative method is 
to use the one- dimensional velocity dispersion of each massive halo as the rms for random 
numbers. 

We compute ui at 200 kpc comoving; ai is very ffat at these scales (see Paper I). We 

label this quantity as a^^^; MH is used to denote the original massive halo. We label 
the i^^ (for i — x, y, z) component of the center-of-momentum velocity of the massive 

halo as vf^^\ We then compute three gaussian random numbers, r^, with mean zero 

and a one-dimensional standard deviation a^^^ for each cluster member. We define the 
velocity of the added cluster member as 

[cluster member] = v^^^^ + 0^^'^ri , (6.1) 
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for some constant P < 1 discussed next. 

The quantity (3 is the ratio of "galaxy temperature" to the virial or gas temperature 
(see Sarazin 1988; Evrard 1990). The "galaxy temperature" is a measure of the kinetic 
energy of the galaxies and the gas temperature is directly related, in hydrostatic equilib- 
rium, to the gravitational potential well. Observational estimates yield (3 ~ 0.8 (Evrard 
1990) with a range 0.4 to 1 (Table 2 of Sarazin 1988). 

We show (jp in Figure 14a and 14b using this method to assign velocities to the added 
cluster members with /? = 1, 0.8, and 0.25. We use the mass-to-light method to break 
up the halos with Vdrc > 350 km s-^ We use M/C = 250 in Figure 14a, M/C = 500 in 
Figure 14b, and in both cases we consider halos with V^irc > 150kms~^; these values, 
and P = 0.25, are chosen specifically to give low estimates of Op. We want to know 
how much we need to "push" the parameters to match the DP estimates of the pairwise 
velocity dispersions. Admittedly, (3 = 0.25 is far below the lowest observational estimates 
(~ 0.4 at best), and is only shown as a final, futile attempt to save CDM! The solid 
curves are for the mass. The /3 = 1 cases are comparable to the M/C = 250 cases using 
the random sampling method of Figure 13. Here they are slightly lower because we show 
halos with Vdrc > 150kms~^ rather than for V^circ > 250kms~^ used in Figure 13. We 
conclude that even small (3 cannot save erg ^ 0.7. The case erg = 0.5 still has pairwise 
velocity dispersions that are high compared with the observations for M/C — 250 and 
the model requires /3 ^ 0.25 which is extremely small compared with observed estimates. 
The same conclusion holds for M/C — 500 at dg = 0.5 except that the (3 — 0.25 case is a 
reasonable match to the observed pairwise velocity dispersions. However, as mentioned 
earlier, /? = 0.25 is far below any observed estimate from real clusters. 

DP computed pairwise velocity dispersions, crp(r), with and without the removal 
of three clusters: Virgo, Coma, and A1367. The triangles in Figs. 14abc are the DP 
numbers computed without these three clusters, yet their effect is small except for the 10 
Mpc bin where the numbers are unreliable. However, the effect is not small for ag = 1 
CDM as we now show. The CDM model is plagued with far too many massive halos (see 
Paper I) and too many rich clusters (see § 5 of this paper). In CDM16, we find 17 halos 
with V^irc > 1000 kms~^ at Ug = 1 (involving 395 galaxies with Vcirc ^ 150kms~^ using 
the mass-to-light method with M/C = 500). We find 2 halos with Vcirc > 1000 km s~^ 
at (jg = 0.7 (involving 57 galaxies with Vcirc > 150kms~^, again using the mass-to-light 
method with M/C = 500). Last, we find no halos with Vcirc > 1000 kms"-*^ at erg = 0.5. 

We now compute (Tp{r) without the inclusion of halos (or added members) with 
V'circ > 1000 kms"-*^. The results are shown in Fig. 14c for halos with Vcirc > 150 kms"-*^, 
P = 0.8, and M/C = 500 applied to halos with Vcirc > 350 km s"^ (cf. Fig. 14b which 
includes all clusters). The effect is substantial at erg = 1; nevertheless, crp(r) is still too 
large by a factor ~ 2 compared with observations for r ^ 1 Mpc. In order to significantly 
reduce crp(r) at ug = 1 we must 1) assume a ridiculously small P and 2) remove an 
extreme number of rich clusters. Even so, these effects are not strong enough to reconcile a 
erg >, 0.7 CDM universe with observed small-scale pairwise velocity dispersions of galaxies. 

As a final comment, we note some recent work that is relevant to the formation 
of galaxies in massive clusters in the CDM model. Katz & White (1993) (see also 
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Evrard et al. 1994) have performed a gas dynamical CDM simulation with = 1, 

Hq = 50kms~^ Mpc~^, and erg = 0.4. They simulated a volume of space containing a 
massive halo found from a previous dark matter only simulation. The object at z = 
(2; = 1/a — 1; a = 1 at a"8 = 0.4) had a mass of 1.83 x 10^^ M0 and a circular velocity of 
945kms~^; it formed from the merging of two massive subclumps. The gas dynamical 
simulation, assuming a 10 to 1 ratio of dark matter mass to gas mass, was evolved to 
z = 0.13. During the course of the simulation eight galaxies formed, but hy z = 0.13 only 
four galaxies had survived the merging process. (Each of these four galaxies had a cold 
gas mass exceeding 1.9 x 10^^ M©.) Estimates of (3, with these limited statistics, were 
/? ~ 1 at z = 0.6 and /? ~ 0.4 at z = 0.13. 

This work is interesting because it demonstrates that some galaxies can survive the 
merging process. However, this fact does not solve the problems demonstrated in this 
paper. Our mass-to-light method predicts that a 1.83 x 10^^ Mq object should contain 
4 halos with V^irc > 250 km s"^ if M/C ~ 500 = lOOO/i. This is a factor of 5 too few 
compared with more typical M/C ^ 200/i; see Trimble (1987). 

Although the gas dynamical simulation of Katz & White (1993) demonstrated that 
some galaxies can survive the merging process in a single massive halo, it did not demon- 
strate that the Q ~ 1 CDM model can successfully make clusters of galaxies with rea- 
sonable mass-to-light ratios. Furthermore, it did not demonstrate that CDM with gas 
dynamics can solve the problems we have found in this paper. We leave open the possi- 
bility that a full scale CDM simulation with gas dynamics might significantly alter the 
distribution of luminous galaxies compared with dark matter halos. However, it is dif- 
ficult to imagine how this could avoid the problems associated with having too many 
massive halos where galaxies are sure to form as seen in the cosmological gas simulations 
of Katz et al. (1992). 

On the other hand, we have found that M/C = 500 at ag = 0.5 might solve some 
of the problems with the models. The numbers of halos and group properties were in 
good agreement with the observations. However, the correlation length (ro ~ 6 Mpc) fell 
short of the observed value ro = 10 Mpc. We found in this section that the velocities 
for M/C ^ 250 at as = 0.5 are marginally consistent with the observed pairwise velocity 
dispersions and are in good agreement with the observed pairwise velocity dispersions for 
/3 ;S 0.25 and M/C — 500. If CDM is to survive on small scales, nature must conspire to 
hide a lot of dark matter. 
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7. CONCLUSIONS 



There appears to be no linear normalization of the power spectrum for the f2 = 1 
CDM model that can simultaneously match the observed numbers, the spatial clustering, 
and the pairwise velocity dispersion of resolved dark matter halos. The problems are 
especially serious for the large amplitude (us ~ 1.0) implied by the recent COBE-DMR 
anisotropy results. 

We must break up the massive halos if our catalogs are to contain groups of galaxies 
like the observed universe. If we study the models without breaking up the massive halos, 
then we find that the two-point correlation function turns over on small scales and the 
correlation length is too small except for erg ~ 1 where the turnover on small scales is 
particularly severe. We also find that the pairwise velocity dispersions constrain < 0.5 
despite the fact that there is a velocity bias of a factor ~ 2. 

We paid considerable attention to massive halos which might represent groups of 
galaxies. Breaking up these massive halos into groups of galaxies removes the turnover of 
the two-point correlation function on small scales and it increases the correlation length 
on larger scales. Unfortunately, the groups do more harm than good unless we assume 
very high mass-to-light ratios. They significantly increase the number of halos, they 
give the wrong shape of the two-point correlation function, they significantly increase 
the pairwise velocity dispersions, and they make groups that are too rich for reasonable 
mass-to-light ratios. Our estimates constrain the models to very high mass-to- light ratios 
~ lOOO/i, although the precise values are uncertain. Factors such as 1) j3, 2) how much of 
the bound mass should be used to estimate the group luminosity (i.e. the mass out to a 
given distance from the group center), and 3) variable mass-to-light ratios, all complicate 
the interpretation of our estimated M/C The combined uncertainty can be as much as 
a factor of ~ 2. However, there is increasing evidence from X-ray studies of clusters that 
dark matter is not hidden in the outskirts of galaxy clusters (e.g. Sciama, Salucci, & 
Persic 1992 and references therein). 

The problems associated with ag ~ 0.4 are clear. In agreement with White et al. 
(1987) we found that we needed to restore halos in massive systems to get the required 
two-point correlation length for as — 0.4. However, the fact that the model then had a 
factor ~ 3 too many halos and produced the wrong shape of the two-point correlation 
function is a serious shortcoming of the model. We also studied models with erg > 0.5. 
In agreement with Couchman & Carlberg (1992) we found a velocity bias of a factor 
~ 2 for (Tg = 1. However, restoring the merged halos in massive systems which have 
high velocity dispersions significantly increased the pairwise velocity dispersions. We can 
rule out CTg ^ 0.7; even erg = 0.5 required a ratio of galaxy to virial temperature P ^ 0.25 
which is too small compared with observed estimates. Removal of the most massive halos 
(with Vcirc > 1000 km s~^) can reduce pairwise velocity dispersions, but the effect is too 
little to save CDM with erg = 1. 

If we live in an O = 1 universe, nature (or clever humans) must learn to hide 
large amounts of dark matter. Gas dynamical simulations probably will not solve the 
problems we have found unless our assumptions regarding sites of galaxy formation and 
galaxy luminosities from the dark matter are significantly wrong. 
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FIGURE CAPTIONS: 



FIG. 1: a) Linear theory predictions for erg as a function of box size Amax- The power 
spectrum (Holtzman 1989, 5% baryons) is normahzed so that as — 1 when Amax — oo. 
b) Linear theory predictions for three-dimensional pairwise velocity dispersions, Ui,, as a 
function of box size, Amax, and separation r, for four values of Amax- 

FIG. 2: log^Q ^(r) (mass) versus log^^g where r is measured in comoving Mpc for various 
simulations in boxes ranging from 51.2 Mpc on a side to 400 Mpc on a side. The five 
128^ particle PM simulations (-R1/2 = 280 kpc, 51.2 Mpc box) are averaged together 
(solid curves) with la error bars. The other simulations are CDM11(128'^,102.4,560; dot- 
dashed curves), CDM16(1443, 100,85; short-dashed curves), and CDM10(1283,400,2188; 
long-dashed curves). Note that ^(r) is underestimated in the 51.2 Mpc boxes. 

FIG. 3: (Tp(r) for the mass for the cases considered in Figure 2. 

FIG. 4: Average log]^Q^(r) for resolved halos with a) T^irc > lOOkms"^, b) Vdrc > 
200kms~^, and c) Vdrc ^ 250kms~^ from two 64"^ particle, P'^M simulations with e =40 
kpc and 51.2 Mpc boxes. The solid curves are for the mass {la error bars are shown 
in the bottom panel) and the solid lines are the observed ^. The results are shown for 
halos found with EOF (/=0.1; short-dashed curves) and FOF {1=0.2; long-dashed curves) 
and DENMAX (512^ grid; dot-dashed curves). In the bottom panel we also show points 
(with la error bars) for catalogs with the massive halos broken up. 

FIG. 5: (Tp(r) at as = 0.7 from CDM12(643,51.2,52). The solid curve is for the mass. The 
dot-dashed curve is for the mass with the particles from the two largest halos removed. 
The result for the DENMAX halos with Vdrc > 192kms~-^ is the long-dashed curve. 
The result using the velocity of the maximally bound particle in a halo rather than the 
center-of-momentum of the halo is the short-dashed curve. 

FIG. 6: logio^(r) for CDM16(1443, 100,85) where the "galaxies" are tagged from the 
initial density field smoothed with a gaussian smoothing radius of Rs = 550 kpc comoving 
{u = 2.6 dot-long-dashed curves; u = 3.0 short-dashed curves) and 880 kpc comoving 
{u = 2.0 long-dashed curves; v = 2.5 dot-short-dashed curves). The curved solid curve is 
for the mass and the solid line is the observed ^. 

FIG. 7: A massive halo at erg = 0.5 from CDM16. a) The massive halo [with mass 2.1 x 
10^^ Mq and V;irc(200 kpc) = 666 km s"^] found by DENMAX; note that DENMAX fails 
to reveal some visually distinct substructure, b) Peak particles that are bound members 
of the massive halo. Circles are for Vdrc > 250kms~^ and squares are for 200kms~^ < 
Vciic < 250 kms"-*^. c) Halos that formed at erg = 0.2 that fell into the massive halo. 
Circles are for V^irc > 250kms~^. d) Halos added to the massive halo assuming a mass- 
to-light ratio of 125; positions are randomly sampled from the massive halo; circles are 
for Vdrc > 250kms~^ and squares are for 192kms~^ < Vdrc < 250kms~^. 

FIG. 8: logio^(r) for DENMAX halos from CDM16. The sohd line is the observed ^. 
The curved line is for the mass. Here there is no special treatment of massive halos. 
The results are shown for Voire > lOOkms"^ (dot-long-dashed curves), Vdrc > 150kms~^ 
(short-dashed curves), V^irc > 192kms~^ (long-dashed curves), and Vdrc > 250kms~^ 
(dot-short-dashed curves) . 
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FIG. 9: logio C(^) for DENMAX halos from CDM16 at as = 0.5. We broke up massive 
halos with Vdrc ^ 350kms~^ using halos that feU in from ag = 0.2 (dot-long-dashed 
curve) and erg = 0.3 (short-dashed curve). We also show results assuming a mass-to-light 
ratio of 125 where we randomly sample the positions (long-dashed curve) and where we 
put all added halos on top of each other (dot-short-dashed curve). The solid line is the 
observed ^. 

FIG. 10: logio ?(^) for DENMAX halos with Kirc > 250 kms"^ from CDM16. We break 
up the massive halos {Vdrc > 350kms~^) using mass-to-light ratios: 50 (dot-long-dashed 
curves), 125 (short-dashed curves), and 250 (long-dashed curves). The results without 
break-up are shown as dot-short-dashed curves. The solid curves are for the mass and 
the solid lines are the observed ^. 

FIG. 11: Cumulative number of galaxies in groups with three or more members; see 
eq. (5.1). We break up the massive halos using the mass-to-light method with M/C — 125 
(top lines), 250, and 500; the higher curves are for the smaller values of M/C. The lowest 
curves are results without breaking up the massive halos. We assume dp/p = 80 to 
determine the FOF linking length. The results are for a 6° wedge out to R = 240 Mpc 
with a magnitude limit B{0) — 15.5. (Note well, we divide the numbers from a 12° wedge 
in the simulation by two. This explains why the jumps in the solid histograms are half 
the value implied by A^mem-) The RGH results for a 6° wedge are shown for comparison 
as dashed histograms. 

FIG. 12: cTp for DENMAX halos from CDM16. Here there is no special treatment of 
massive halos. We show circular velocity cut-offs of lOOkms"^ (dot-long-dashed curves), 
150kms~^ (short-dashed curves), 192kms~^ (long-dashed lines), and 250 kms~^ (dot- 
short-dashed curves). The observed estimates are shown as open symbols for various 
modeling parameters from Davis & Peebles (1983). The solid curves are for the mass. 

FIG. 13: CTp for DENMAX halos with Kirc > 250 kms"^ from CDM16. We break up the 
massive halos (Kirc > 350 kms"-*^) using the mass-to-light method: M/C = 50 (dot-long- 
dashed curves), M/C = 125 (short-dashed curves), M/C = 250 (long-dashed curves). 
(The positions and velocities of the added halos are assigned using random sampling.) 
The dot-short-dashed curves are without break-up. The solid curves are for the mass. 

FIG. 14: (Tp for DENMAX halos with VcWc > 150 km s"^ from CDM16. The massive 
halos (Fcirc > 350kms-i) are broken up with a) M/C = 250 and b) M/C = 500. 
The solid curves are for the mass. The results without break-up are the dot-short- 
dashed curves. The velocities of added members are generated from the central velocity 
dispersions and /? (eq. (6.1)). The M/C curves are for /3 = 1.0 (dot-long-dashed curves), 
P = 0.8 (short-dashed curves), and /? = 0.25 (long-dashed curves). c)(7p(r) for CDM16 
halos with Vdrc > 150kms~^, removing those halos (and their progeny) with V^circ > 
1000 kms~^. Halos with 350kms~^ < Vdrc < 1000 kms~^ are broken up with P — 0.8 
and M/C = 500. 
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